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We study the classical geodesic motions of nonzero rest mass test particles and photons in 
five-dimensional warped product spaces. We show that it is possible to obtain a general picture 
of these motions, using the natural decoupling that occurs in such spaces between the motions 
in the fifth dimension and the motion in the hypersurfaces. This splitting allows the use of 
phase space analysis in order to investigate the possible confinement of particles and photons to 
hypersurfaces in five-dimensional warped product spaces. Using such analysis, we find a novel 
form of quasi-confinement which is oscillatory and neutrally stable. We also find that this class 
of warped product spaces locally satisfy the Z2 symmetry by default. The importance of such 
a confinement is that it is purely due to the classical gravitational effects, without requiring the 
presence of brane-type confinement mechanisms. 
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I. INTRODUCTION 

Motivated by string inspired theories [H, a great 
deal of effort has, over the recent years, gone into the 
study of higher dimensional general relativity (GR) the- 
ories. Examples of such theories include noncompacti- 
fied (n + l)-dimensional GR [|| (where the observable 
universe is assumed to lie on an embedded 3 + 1 (4D) 
hypersurface), Kaluza-Klein theory [i| (where the extra 
dimensions are compactified) and the braneworld models 
inspired by String/M-theory @, d, @. 

A central issue in such theories has been how to explain 
the fact that the observable universe is (very nearly) con- 
fined to a particular hypersurface. In most braneworld 
models, stringy effects are invoked to argue that in low 
energy regimes particles are restricted to a special 3 + 1 
brane hypersurface, which is embedded in a higher (usu- 
ally five)-dimensional bulk, while the gravitational field 
is free to propagate in the bulk |l| . 

In the case of other higher-dimensional GR theories, 
however, such confinement needs to be explained. This 
has provided strong motivation for the consideration of 
the so-called warped product spaces [7j and the investi- 
gation of their geometrical properties. We recall that the 
geometrical considerations of warped product spacetimes 
predate their applications in these higher-dimensional 
contexts and were first considered by Bishop and O 'Neill 
[H in the context of 3 + 1 general relativity. 

Geodesic motion in connection with higher- 
dimensional theories has been considered by a number 
of authors in recent years [H, [l(| EH- Our interest in 
studying such motions specifically in the context of 



warped product spaces arises from the realization of the 
fact that there is a strong connection between the 5D 
geodesies and the analytical form of the warping func- 
tion. Indeed, as we shall see in this paper, it is possible 
to make a general qualitative analysis of the behaviour 
of massive particles and photons in the fifth dimension 
from the knowledge of the warping function. This can 
in turn allow questions concerning the possibility of 
hypersurface confinement as well as the stability of the 
confined motions with respect to small perturbations to 
be treated with great generality. 

An assumption usually made in the braneworld scenar- 
ios is that the branes possess an infinitesimal thickness. 
In this sense one may talk about 6 -confinement of mas- 
sive particles and light rays. There are also smooth ver- 
sions of this scenario which consider the possibility that 
branes may have a finite thickness [H, [H, EH • In that 
case one may talk about quasi-confinement. 

Here we shall demonstrate a novel example of the lat- 
ter in the context of 4 + 1-dimcnsional warped product 
spaces which we shall refer to as oscillatory confinement, 
whereby the particles in the four-dimensional hypersur- 
face can oscillate about the hypersurface while remain- 
ing close to it. As we shall show below such behaviour 
can indeed occur for large classes of bulks which possess 
warped product geometries. Our analysis is fairly general 
and the results may be applied to different examples of 
warped product spacetimes considered in the literature. 

The paper is organized as follows. In Section 2 we write 
the geodesic equations for warped product spaces and 
consider the parts due to 4D and the fifth dimension sep- 
arately. We then show that the equation that describes 
the motion in 5D decouples from the rest. We proceed 



2 



in Section 3 to rewrite the geodesic equation in the fifth 
dimension as an autonomous planar dynamical system. 
We then employ phase plane analysis to study the motion 
of particles with nonzero rest mass (timelike geodesies) 
and photons (null geodesies) respectively. Such qualita- 
tive analysis allows very general conclusions to be drawn 
about the possible existence of confined motions and their 
stability in the neighbourhood of hypersurfaces. We con- 
clude with a discussion of further applications in Section 
4. 

Throughout Latin indices take values in the range 
(0,1, ...,4) while Greek indices run over (0,1,2,3). 



II. WARPED PRODUCT SPACES AND 
FIVE-DIMENSIONAL GEODESICS 

Warped geometries are at the core of the Randall- 
Sundrum models. Technically we can define a warped 
product space in the following way. Let (M m ,h) and 
(M n ,k) be two Riemannian (or pseudo-Riemannian) 
manifolds of dimension m and n, with metrics h and 
k, respectively. Suppose we are given a smooth function 
/ : M™ — > R (which is referred to as the warping func- 
tion) . Then we can construct a new Riemannian (pseudo- 
Riemannian) manifold by setting M = M n x M m and 
defining a metric g = e 2 ^h®k. In this paper we shall take 
M = M 4 x K, where M 4 is a four-dimensional Lorentzian 

manifold with signature (H ) (referred to as the 

(3+l)-dimensional spacetime) . In local coordinates {y a } 
the line element corresponding to this metric will be de- 
noted by 

dS 2 = 9abdy a dy b . 

Let us now consider the equations of geodesies in the 
five-dimensional space 
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dy b dy c 
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where A is an affine parameter and < - 5 - ) r a c are the 5D 
Christoffel symbols of the second kind defined by ^'T^. = 
\g ad {9db,c + gdc,b - 9bc,d)- Denoting the fifth coordinate 
y 4 by I and the remaining coordinates y^ (the "space- 
time" coordinates) by i.e. y a = (x M , I), it is not diffi- 
cult to show that the " 4D part" of the geodesic equations 
(fTj) can be rewritten in the form 



fx^ (4) dx^dx^_ 
d\ 2 + ^ dX dX 9 ' 



(2) 



and = \g» v (g va ,p + g v p, a ~ g a p,v)- We shall as- 

sume that the five-dimensional manifold M can be foli- 
ated by a family of hypersurfaces S defined by the equa- 
tion / = constant. Then the geometry of each hyper- 
surface, say I = lo, will be determined by the induced 
metric 

ds 2 = g a p{x, lo)dx a dx °. 

Therefore the quantities which appear on the left- 

hand side of Eq. © may be identified with the Christoffel 
symbols associated with the induced metric in the leaves 
of the foliation defined above. 

We shall consider the class of warped geometries given 
by the following line element 



dS 2 = e 2f h a0 dx a dx 13 - dl 2 , 



(4) 



where / 
easy to see [2C 



f(l) and hap — h a p{x). For this metric it is 
at C 6 >I* = and = \g^g^ vA = 

f'Sft, where a prime denotes a derivative with respect to 
I. Thus in the case of the warped product space (UJ) the 
right-hand side of Eq. © reduces to ^ = -2/'^^ 
and the 4D part of the geodesic equations becomes 



d 2 x^ 
~dX 2 ~ 



_(4) p 



dx a dx? 
al3 ~dX~dX 
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dx^ dl 
~dX~dX' 
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Likewise the geodesic equation for the fifth coordinate I 
in the warped product space becomes 



d 2 l , 2 - dx a dx 



(6) 



By restricting ourselves to 5D timelike geodesies 
1 ) we can readily decouple the above equa- 



tion from the 4D spacetime coordinates to obtain 



= 0. 




(7) 



Similarly, to study the motion of photons in 5D, we must 
consider the null geodesies {^)ab^§x^d\ = m which 
case Eq. © becomes 



(8) 



where 



_(5)pp 
1 44 



g 11 (g4a,f3 + g^ji.a ~ ffa/9,4) 
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_,, dx a dl 
ai dX dX 



dx a dx p 
~dX~dX' 



(3) 



Equations and © are ordinary differential equa- 
tions of second-order which, in principle, can be solved if 
the function /' = f'(l) is given. Conversely, if I = l(X) is 
known then / = /(/) can be determined (up to a constant 
of integration) provided that we can write A = X(l). 
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III. QUALITATIVE ANALYSIS OF MOTION IN 
THE FIFTH DIMENSION 

We recall that the hypersurface sets the boundary con- 
ditions for the warping function. If the warping function 
/ = f(l) is given a priori, then we can obtain a great deal 
of information about the motion in the fifth dimension 
by performing a qualitative analysis of Eqs. {7J and l[5]). 
without the need to solve these equations analytically. 
To do this we define the variable q = 4r and consider 
instead of ([7]) or the autonomous dynamical system 



dl 

dX = q 

fx = < 9 > 

with P(q, I) = — f'(e + q 2 ), where e = 1 in the case of ([7]) 
(corresponding to the motion of particles with nonzero 
rest mass) and e = in the case of ([8]) (corresponding 
to the motion of photons). In the investigation of dy- 
namical systems a crucial role is played by their equilib- 
rium points, which in the case of system ^ are given 
by j[ = = . The knowledge of these points to- 
gether with their stability properties allows a great deal 
of information to be gained regarding the types of be- 
haviour allowed by the system. Thus if for a particular 
geodesic we know t as a function of the parameter A, i.e. 
t = t(X), then since ^ ^ 0, knowing the behaviour of 
I = l(X) from qualitative analysis will enable us to de- 
duce the time evolution of the motion of a particle (or 
light ray) in the fifth dimension. 

A. The 5D motion of particles with nonzero rest 
mass near the hypersurface 

We begin by considering the case of nonzero rest mass 
particles, whose motion in the fifth dimension is governed 
by the dynamical system 

dl 

Tx = q 

(10) 

The equilibrium points of (fT0|) are given by q — and the 
zeros of the function f'(l) (if they exist) which we gener- 
ically denote by ^o- These solutions, pictured as fixed 
points in the phase plane, correspond to curves which 
lie entirely in a hypersurface £ of our foliation (since 
they have I = constant). It turns out that these curves 
are timelike geodesies with respect to the hypersurface 
induced geometry. Indeed, the existence of equilibrium 
points of the dynamical system ([9]) has the important 
consequence that the geodesies of the bulk and that of the 
hypersurface (Iq = constant) coincide. This can be seen 
directly from the equations © and ©, and is perfectly 



consistent with a well-known theorem of differential ge- 
ometry according to which the geodesies of a Riemannian 
space and that of a particular hypersurface coincide (i.e. 
confinement occurs) if and only if the extrinsic curvature 
of that hypersurface vanishes (see e.g. [l6| where such 
hypersurfaces are referred to as totally geodesic). Note 
that for the metric ^ the extrinsic curvature of the hy- 
persurfaces £ is given by Sl a p = —f'e 2 fh a p(x), which 
is clearly zero at the equilibrium points, where /' = 0. 
Thus in the absence of equilibrium points (i.e. /' ^ 0) 
the extrinsic curvature does not vanish and according to 
the above theorem the bulk and hypersurface geodesies 
cannot be identical. 

To obtain information about the possible modes of 
behaviour of particles and light rays in such hypersur- 
faces, it is important to study the nature and stability 
of the corresponding equilibrium points. This can be 
done by linearising equations (|10p and studying the 
eigenvalues of the corresponding Jacobian matrix about 
the equilibrium points. Assuming that the function /'(/) 
vanishes, at least at one point lo, it can readily be shown 
that the corresponding eigenvalues are determined by the 
sign of the second derivative f"(lo), at the equilibrium 
point, and the following possibilities can arise for the 
equilibrium points of the dynamical system (|10[) : 

Case I. If f"(lo) > 0, then the equilibrium point 
(q = 0, 1 = Iq) is a center. This represents a case in 
which the solutions near the equilibrium point have the 
topology of a circle. In that case the phase portrait 
consists of closed curves describing the motions of 
particles oscillating about the hypersurface £ (I = la) 
indefinitely (see the top panel in Fig. [TJ) . The amplitudes 
of the oscillations will depend on initial conditions. 
We note that the existence of such cyclic motions is 
independent of the ordinary 4D spacetime dimensions, 
and, except for the conditions f'{lo) = and f"(lo) > 0, 
the warping function f(l) remains completely arbitrary. 
The presence of such center equilibrium points with 
the consequent quasi-confinement of particles can be 
viewed as providing examples of almost totally geodesic 
hypersurfaces (see the lower panel in Fig. [I}. 

Case II. If f"{lo) < 0, then the point (q = 0, 1 = lo) is 
a saddle point. In this case the solution corresponding to 
the equilibrium point E is highly unstable and the small- 
est perturbations will lead to exponential divergence of 
the solutions and hence the particles - from the hyper- 
surface (See Fig. 2). The only (measure zero) exception 
is provided by the stable manifold of the point E, repre- 
sented by the lines AE and BE, along which particles are 
attracted towards the brane. An example of this highly 
unstable "confinement" is provided by the warping func- 
tion [m 

f(l) = -61ncosh(d), (11) 

where b and c are positive constants. Here we have 
a unique equilibrium point at I = and it is easy to 
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FIG. 1: Depicted in the top panel is the center equilibrium 
point E in the case where f"(la) > 0. In this case the particles 
oscillate about the hypersurface I = Iq. This provides a quasi- 
confinement mechanism whereby massive particles enter and 
leave the hypersurface E indefinitely (lower panel). 






FIG. 2: When f"(lo) < the equilibrium point E is a saddle. 
In this case confinement is highly unstable. 



out in ref. [ll[ this is a highly unstable confinement in 
the sense that any small transversal perturbations in the 
motion of massive particles along the brane will cause 
them to be expelled into the extra dimension. However, 
this case lies somewhat outside the scope of the metrics 
that we are considering here as in this case the warping 
function is not smooth (the first derivative of f(l) with 
respect to the extra coordinate is not continuous at 
1 = 0). In the braneworld scenario, as is well known, 
stable confinement of matter fields is possible at the 
quantum level if we take into account interaction with a 
scalar field 17]. In a purely classical picture, however, 



one would require mechanisms other than geodesic 
confinement in order to constrain massive particles to 
move on hypersurfaces in a stable way. 



verify that /"(0) < in this case. In this connection we 
note that for large values of I the warping function (fTT]) 
approaches that of the Randall-Sundrum metric 0]) 



Case III. If f"(lo) — 0, then both eigenvalues are 
zero, which corresponds to a degenerate case. To 
carry out a qualitative analysis of the solutions near 
the equilibrium point, in this case, would require the 
knowledge of the third (or higher) derivative of the 
function f(l). We shall not consider this case in its 
generality. However, if f(l) is constant (case IV), then 
we can easily draw a global phase portrait of the system. 



ds 2 = e- 2k ^r laP dx a dx !3 - dl 2 , 

where k is a constant. In this case f'(l) = according 
to whether I is positive or negative. Hence for / =/= 
there exist no equilibrium points, and therefore no 
confinement of particles purely due to gravitational 
effects. However, in this thin wall limit if we take the 
extrinsic curvature of the brane / = as being zero, then 
(jTUJ) implies confinement. Nevertheless, as was pointed 



Case IV. If f(l)= constant, then /'(/) vanishes for all 
values of I, which implies there is an infinite number of 
non-isolated equilibrium points, that is, we have a line 
of equilibrium points (q = 0). Perturbations along this 
line are neutrally stable, which implies that particles 
placed on any of the hypersurfaces of the foliation I = 
constant, will remain there so long as they do not receive 
a transversal velocity to the hypersurface. 
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Case V. If there are no equilibrium points, that is, the 
warping function f(l) does not have any turning points 
for any value of I, then we cannot have confinement of 
classical particles to hypersurfaces solely due to gravita- 
tional effects. An example of this situation is illustrated 
by the warping function f(l) = | In (A/ 2 /3) considered in 
Ref. [3. 

We note that even when there are no equilibrium 
points, the picture of dynamics can still be obtained if 
f(l) is known, by obtaining a first integral of the system 
©. This can easily be done by writing this system as 
the first-order differential equation 



dl 

dq 



which can be readily integrated to yield 



f(l) = -\ny/e + q 2 + K, 
where K is a constant of integration. 



B. The 5D motion of photons near the 
hypersurface 



(12) 



(i) f(l) is an increasing monotonic function for I > 
which approaches the limit f(l) — > 00 as I — » 00 
(see the top panel of Fig. [3]). This includes the case 
of the warping function considered in ([H[). 

(ii) f(l) is a decreasing monotonic function for I > 
that approaches the limit f(l) — > —00 as I — ► 00, 
which is the case of the warping function given by 
Ijlip (see the lower panel of Fig. |5J). Note the inter- 
esting behaviour of photons that are not confined: 
if we set them in motion towards the hypersurface 
E(Z = 0) then after reaching a minimal distance 
from E they bounce back to infinity, showing clearly 
that they are repelled by E; 

(iii) f(l) is constant, in which case we have the same 
phase diagram as in the case of massive particles. 

Finally it is worth noting that in the case of null 
geodesies, the equilibrium points of the Eqs. (fTB")) do 
not require /' = 0, which is a consequence of the fact 
that the above mentioned theorem [l(| does not apply to 
null geodesies. 



IV. CONCLUSIONS 



Finally we examine the motion of photons (null 
geodesies) and ask whether they can be confined to or 
a neighbourhood of a hypersurface. In this case the dy- 
namical system (|9]) becomes 



dl 

dX 
dq 

dX 



= q 



= -/v 



(13) 



The equilibrium points now are given by q — 0, which 
consists of a line of equilibrium points along the Z-axis, 
with eigenvalues both equal to zero. As a result there 
exist 5D null geodesies in any hypersurface I — constant. 
This demonstrates that confinement of photons to hyper- 
surfaces does not depend upon the warping factor. There 
are also regions of stability and instability with respect 
to small perturbations along the Z-axis. As was pointed 
out in the previous section we can easily obtain a first 
integral of the system (fl3|) which is given by ( IT2|) . In 
the case of the motion of photons (e = 0) this gives 

q = Ae-M 

where A is a constant of integration. Therefore we can 
obtain a global picture of the solutions (fT3|) provided that 
we have some qualitative knowledge of the function f(l). 

Many different cases can arise, depending upon the 
nature of f(l). Here for definiteness and simplicity we 
impose Zi symmetry on the geometry of the bulk. The 
following are some possible examples: 



We have studied some aspects of the motion of massive 
particles and photons in five-dimensional warped prod- 
uct spaces. Spaces of this type have received a great deal 
of attention over the recent years mainly in connection 
with the braneworld scenarios. Our treatment has been 
geometrical and classical in nature and we have consid- 
ered confinement to hypersurfaces purely due to classical 
gravitational effects. 

Employing the splitting that naturally occurs in such 
spaces between the motion in the fifth dimension and the 
remaining dimensions, we have found, using phase plane 
analysis, a novel form of quasi- confinement which is neu- 
trally stable. In a sense this amounts to a generalisation 
of the S -confinement of particles in thin branes. An im- 
portant reason to consider models outside the framework 
of thin branes is that there is a characteristic minimum 
length given by the string scale [l(| [H, EE 03 • 

The importance of the type of confinement discussed 
here is that it is due purely to the classical gravitational 
effects, without requiring the presence of brane type con- 
finement mechanisms. 

As an example, we have found that in the case of 
the Randall-Sundrum metric gravity alone is not capable 
of confining massive particles, although photons may be 
confined. As is well known, the particle confinement in 
this scenario is achieved by means of a scalar field 0]. 

An assumption usually made in braneworld scenarios, 
which plays an important role in connection with con- 
finement in these models, is that of Z2 symmetry. It is 
therefore important to take a closer look at the models 
we have found, which are capable of providing successful 
confinement, from this perspective. 
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FIG. 3: Depicted in the top panel is the phase portrait rep- 
resenting the motion of photons in presence of Z2 symmetry 
with an increasing monotonic function f(l) for / > and 
/(/) — > 00 as / — > 00. There is confinement of photons at 
any hypersurface E of the bulk foliation I — constant. The 
bottom panel shows the corresponding phase portrait for the 
cases where f(l) is a decreasing monotonic function for I > 
that approaches the limit f(l) — » —00 as I — » 00. Photons that 
are not confined are repelled by the hypersurface E (I — 0). 



Recalling that Z2 symmetry implies 
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it is easy to check that this condition is satisfied in cases 
where the hypersurfaces are located at the turning points 
I = Iq of the functions /. This follows from the fact that 
in such cases /', and hence VL a p = — f'e 2 * h a p(x), have 
different signs on cither sides of the hypersurface and 
thus Z2 symmetry is satisfied by default. Importantly 
this includes all the equilibrium points we have found in 
the massive particle case. 



Finally it is interesting to compare our results with the 
classical theorem, according to which the geodesies of a 
Riemannian space coincide with the geodesies of an em- 
bedded (totally geodesic) hypersurface if and only if the 
extrinsic curvature of that hypersurface vanishes. The 
analysis involved in proving this theorem (l6| . though 
general, does not address the case of null geodesies, as 
it confines itself to positive definite metrics. Nor does 
it provide any information regarding the stability of the 
this coincidence. Our analysis on the other hand, al- 
though restricted to the case of warped product spaces, 
provides information regarding the stability of confine- 
ment to totally geodesic hypersurfaces, through the use 
of dynamical system analysis. 
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